Order Parameter and Scaling Fields in Self-Organized Criticality 



o 

Q 



J2 

o 



CO 

o3 



i 

C 

O 
o 



> 

00 

in 

c3 



-a 
c 

o 
o 



X 
S3 



Alessandro Vespignani 1 and Stefano Zapperi 2 
1 Instituut-Lorentz, University of Leiden, P.O. Box 9506 2300 RA, Leiden, The Netherlands 
2 Center for Polymer Studies and Department of Physics, Boston University, Boston, MA 02215 

(February 1, 2008) 

We present a unified dynamical mean-field theory for stochastic self-organized critical models. 
We use a single site approximation and we include the details of different models by using effective 
parameters and constraints. We identify the order parameter and the relevant scaling fields in order 
to describe the critical behavior in terms of usual concepts of non equilibrium lattice models with 
steady-states. We point out the inconsistencies of previous mean-field approaches, which lead to 
different predictions. Numerical simulations confirm the validity of our results beyond mean-field 
theory. 
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The origin of scaling in nature jy has become in the 
last years a challenging problem in physics. Bak, Tang 
and Wiesenfeld (BTW) Q have proposed self-organized 
criticality (SOC) as a unifying theoretical framework to 
describe a vast class of driven systems that evolve "spon- 
taneously" to a stationary state, characterized by power 
law distributions of dissipation events. Despite the in- 
sights SOC concepts have brought to a number of prob- 
lems, an agreement on the exact definition of SOC has 
still not been reached and the exact meaning of the word 
"spontaneous" is quite unclear. Originally, SOC was as- 
sociated with the absence of tuning parameters, but it 
has been noted |3|,[4| that the driving rate acts as a tun- 
ing parameter in most, if not all, SOC models. This 
ambiguity has hindered the formulation of precise rela- 
tions between SOC and other non equilibrium critical 
phenomena |]||. 

In this letter we reformulate SOC in terms of typical 
concepts of non equilibrium critical phenomena by us- 
ing the dynamical single site mean- field (MF) theory §. 
We provide a general scheme in which the details of dif- 
ferent models are included via effective parameters and 
constraints. We mainly discuss sandpile models j^j with 
and without dissipation |8j , but the formalism can be di- 
rectly applied to other stochastic SOC models, such as 
the forest-fire model |9| . Work is in progress to apply this 
method also to systems driven by an extremal dynamics 
p0| . We find two independent critical parameters, i.e. 
relevant scaling fields, both with critical value equal to 
zero, and just in this double limit criticality is reached. 
We study the behavior of the order parameter and eval- 
uate critical exponents. The results we obtain are in 
contrast with previous MF approaches 0,ff2[. This is 
due to a subtle inconsistency in the way critical param- 
eters have been chosen in previous works. We show that 
some MF exponents are exact also in low dimensional sys- 
tems because of conservation laws. Our prediction are 
confirmed by numerical simulations of two dimensional 
sandpile models. 

Sandpile models are cellular automata with an inte- 
ger (or continuous) variable Zj (energy) defined in a 
d— dimensional lattice. At each time step an energy grain 
is added to a randomly chosen site, until the energy of a 



site reaches a threshold z c . When this happens the site 
relaxes (zi — > Zi — z c ) and energy is transferred to the 
nearest neighbors (z.j —* Zj + y 3 ■ ) . For conservative mod- 
els the transferred energy equals the energy lost by the 
relaxing site {J^Vj = z c), a t least on average. Usually, 
the only form of dissipation occurs at the boundary, from 
which energy can leave the system. With these condi- 
tions the system reaches a stationary state characterized 
by avalanches whose sizes s are distributed as a power 
law P(s) ~ HIS 

In order to simplify the description of these models, we 
can reduce the number of states each site can assume in 
the following way. We divide sites in critical, stable and 
active |14|. Stable sites are those that do not become 
active if energy is added to them. Critical sites become 
active by the addition of energy. Active sites are relax- 
ing and transfer energy, providing an interaction with 
other sites, usually the nearest neighbors (n.n.). In this 
way we have mapped the system in a three states cellu- 
lar automaton (CA) on a d-dimensional lattice |(]|l5| ■ To 
each site i is associated a variable Sj, which can assume 
three different values. A complete set s = {s,} of lattice 
variables specifies a configuration of the system. The 
dynamics is characterized by the operator (s | W \ s ) 
which represents the transition rate from a configuration 
so to a configuration s in a time step t. A well estab- 
lished technique to study these systems is the single site 
mean- field approximation (?]]. Denoting by p a ,Pc and p s 
the average densities of sites in the active, critical and 
stable states respectively, we write the following reaction 
rate equations 



0_ 



Pk = F K {p a , p c , p s ) 



a, c, s. 



(1) 



Because the densities must preserve normalization, two 
of the above equations supplemented with the condition 
Pa + Pc + Ps = lj are enough to describe completely the 
system. The most general way to represent the function 
F K is through the Taylor series of the average densities: 



F K = £ ftp* + £ r/pnPi + 0{pl), (2) 
n n,t 

where the constant term is set to zero in order to get 
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a stationary state. The first order terms are the tran- 
sition rates generated by the external driving fields or 
by spontaneous transitions. The second and higher or- 
der terms characterize transitions due to the interaction 
between different sites. In SOC models, only the active 
state generates a non trivial dynamical evolution, while 
stable or critical sites can change their state only be- 
cause of the external field or the presence of an active 
n.n. site. Since the critical point is identified by p a = 0, 
in correspondence with a vanishing external field, we can 
neglect second order terms in the density of active sites. 
The solutions of the stationary equations (Jj/0 K = 0) are 
function of the effective parameters /™,/™'^, which de- 
pend on the details of the model. It is expected that the 
critical behavior is not affected by the specific values of 
the parameters, while universality classes will depend on 
constraints imposed on the equations, because of symme- 
tries and conservation laws. 

For the seek of clarity, we describe in details the case 
of sandpile models. In this class of systems the only ex- 
ternal field is the flow of energy added to the system. 
We can describe this driving by the probability per unit 
time h that a site will receive a grain of energy. The total 
amount of energy added to the system at each time step 
will be J; m = hL d . The first order terms in F a are the 
transition rates a — * s, c and vice versa, independently of 
nearest neighbor sites: 



J a 



-l 



f s = o • f c 

Ja u > J a 



(3) 



Here we considered that active sites becomes stable with 
unitary rate, stable sites never becomes active and criti- 
cal sites becomes active because of the external field. In 
addition, there is a single interaction term that describes 
the creation of an active site from a critical site due to 
the relaxation of n.n. sites. We can write this term as 
(<7 — e)p c p a , where g is an effective rate that depends on 
the geometry and the energy involved in the relaxation 
process and e is the average energy dissipated in each site 
Jl6[ . We stress that e is present also for fully conserva- 
tive systems, being an effective term due to the boundary 
dissipation jl7|. Considering all these terms we obtain 



Fa 



-Pa + hp c + (g~ e)p c Pa + 0(p 2 a ). 



(4) 



A similar reasoning yields the functions F c and F s . The 
effect of the driving field on stable sites and the inter- 
action between active and stable sites deserves a discus- 
sion. The corresponding terms are contributing to the 
transition rate s — > c. In sandpiles models, the energy 
conservation imposes a local constraint in the rate equa- 
tions. Energy is stored in stable sites until they become 
critical, but only a fraction u of stable sites receiving an 
energy grain contribute to the s — > c process. Therefore 
in this case the reaction rates will be given by the c — > a 
rates multiplied by the factor u. For instance, a stable 
site will receive an energy grain with probability h, but 
only with probability uh it will turn critical. The same 
reasoning holds for the interaction term. 

After imposing stationarity, we get the following dy- 
namical MF equations: 



Pa = hp c + (g - e)p c p a 
p a = uhp s + u(g - e)p s p a 
Pa = 1 - Ps - Pc 



(5) 



where u, g are effective parameters which depend upon 
the particular model, e represents the dissipation and h 
is the driving field. We expect g to be an independent 
parameter of the model, while u has to be obtained self- 
consistently, because it is fixed once the dynamical rules 
of the CA are given. 

After some algebra from Eqs. we obtain a closed 
equation for p a 

u(g - e)p 2 a + (1 + tt(l + h-g+ e))p a - uh = 0. (6) 

We can expand p a (h) for small values of the field h. The 
zero order term in the expansion vanishes and we obtain 
a leading linear term: 



Pa(h) = 



ih 



1 + u — ug + ue 



(7) 



This result has to be consistent with the global conser- 
vation law, which states that the average input energy 
flux Ji n must balance the dissipated flux J ou t- In the 
stationary state the conservation law can be written as 



Jin 



hL° 



J O 



ep a L a 



(8) 



By comparing Eq. [j] with Eq. ^ we obtain that u — 
1/(5—1). In the limit h — > the densities are there- 
fore given by 

Pa = ~, Pc = - + 0(h), p s = g —± + 0{h). (9) 

£9 9 

An estimate of g can be obtained using a random neigh- 
bor approximation, which yields g = 2d for the BTW 
model H or g — 2 for the two level models No- 
ticeably, in the latter case u = 1, i.e. all stable sites are 
sub-critical, as it is expected for a two level model. 

We now discuss the critical behavior of these systems. 
The balance between conservation laws and dissipation 
is essential for the critical behavior of the model, as it 
has been also pointed out in |l£|. The model is critical 
just in the double limit h, e — > 0, h/e — > 0, similarly to 
the forest-fire model || . In analogy with non equilibrium 
phenomena J6|,[l5| , the one particle density of active sites 
is the order parameter and goes to zero at the critical 
point. We can then distinguish several different regimes 
as a function of the parameters. The system has no sta- 
tionary state for h > e, since p a would have to be greater 
than one to satisfy Eq. (||). The model is supercritical 
for h > and e > h, while for h — ► and e > it is sub- 
critical and the dynamics displays avalanches. The phase 
diagram is somehow similar to that of usual second order 
phase transitions, if we replace h by the magnetic field 
and e by the reduced temperature. 

In the supercritical regime the order parameter is lin- 
ear in h 



Pa 



5=1. 



(10) 
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This is analogous to the MF results obtained for contact 
processes and other non equilibrium CA [|6|j7|Jl5[| , but it 
is in contrast with previous MF approaches for sandpilc 
models 0,0, which yielded S = 2 @. This incor- 
rect result is due to an inconsistency present in those 
studies. The scaling is expressed in terms of the average 
energy 9 = piZi which is treated as an independent 
control parameter. As we have just shown, 9 and h are 
not independent. Moreover, 9 can not be considered as 
the control parameter even for h = 0, since it does not 
determine completely the state of the system: the same 
value of 9 describes several states corresponding to dif- 
ferent values of densities pi. This is a typical property of 
CA with multiple absorbing states [HJ . In analogy with 
non equilibrium CA it is possible to define several other 
exponents characterizing the supercritical regime for this 
class of models. The complete set of exponents and their 
scaling relations will be reported elsewhere Jl6| . 

In the subcritical regime , the behavior of the system 
is dominated by the dissipation. This can be seen by 
studying the susceptibility 



X 



_ dp a {h) 
dh 



(11) 



which diverges for e = 0. The system is in a subcritical 
state for any value of e different from zero. The critical 
behavior is thus characterized by the scaling laws \ ~ e~ 7 
and £ ~ e _l/ , where £ is the characteristic length. 

We can use these exponents to characterize the con- 
servative sandpile model, since our MF analysis treats 
both border and bulk dissipation. In conservative sys- 
tems, when the size is increased the effective dissipation 
decreases as e ~ since bulk processes dominate over 
boundary dissipation. At the same time, the character- 
istic length of the avalanches should go like £ ~ L to 
ensure dissipation of energy. This implies the scaling re- 
lation vp, = 1, and that \ ~ L^ 1 . It is also possible to 
show [fLq] , that the susceptibility scales as the average 
avalanche size, and in two dimensions it has been found 
that (s) ~ L 2 exactly for L — > oo [^0|. The same result 
holds also in MF theory if we assume that the dynamics 
is diffusion like Jl|^||. 

Combining all the above results we obtain a first set of 
MF exponents 
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1, fx = 2, v=l/2. 



(12) 



It is worth to remark that it is not possible to define the 
equivalent of an exponent /3, because for h = the order 
parameter is always zero. We emphasizes again that h 
and e are both control parameter responsible for different 
regimes of the model. 

We have derived these exponents using only conserva- 
tion laws, therefore we expect they should describe also 
low dimensional sandpile models. We simulate numer- 
ically the BTW model with finite driving rate h and 
boundary dissipation. We see in Fig [l] that the density 
of critical sites goes to zero linearly with h (S = 1) with 
a slope that increases with the system size as L 2 . This 
is in agreement with the MF theory which predicts that 
the susceptibility scales as L' 17 , with pj = 2. To observe 



more clearly the scaling with dissipation of the sandpile 
model we study the BTW model with periodic boundary 
conditions and fixed dissipation e [^3| . In Fig. || we plot 
the control parameter as a function of h/e. The scal- 
ing predicted by the MF theory (7 = 1) is verified with 
remarkable accuracy, and we note that finite size cor- 
rections are not noticeable, in contrast with the case of 
boundary dissipation p3[ . Finally, the exponent v = 1/2 
has been measured already in a two dimensional dissipa- 
tive sandpile model [||. 

The dynamics in the subcritical regime takes place 
in the form of avalanches. The exponents describing 
avalanche distributions depends on the dimension and 
in general will not agree with the MF results. A com- 
plete characterization of MF avalanche scaling has been 
obtained by using the theory of branching processes [jl8] . 
Here, we reproduce these results in an independent way. 
Following Grassberger and de la Torre H we consider the 
probability that a small perturbation activate s sites (an 
avalanche in the SOC terminology) 



P(s,e) 



(13) 



where s c ~ e 1 / <T is the cutoff in the avalanche size. The 
perturbation decays in the stationary subcritical state as 



Pa (t) ~ f>T{t/t c ). 



(14) 



Here t c denotes the characteristic time which scales as 
t c ~ e~ A . We can obtain these exponents by solving 
Eq. (|l|) for a small perturbation around the subcritical 
stationary state (h = 0) 



dp a , , « 

~gf = ~P a + (9 ~ e )P<*Pv 



(15) 



where, as a first order approximation, we can replace 
p c = l/g. We obtain in this way p a (t) ~ exp(—et/g) 
which implies rj = and A — 1. Introducing the scaling 
laws s c ~ £ D and t c ~ £ 2 , it is possible to derive pl| 
another set of scaling relations 

7=^, D=\ ZV = A, i^ = z; (16) 
a vo va 

from which we get the second set of MF critical exponents 

z = 2, D = 4, r = 3/2, (7 = 1/2; (17) 

in agreement with the theory of branching processes p8| . 
It is worth to remark that the numerical value of these 
exponents is the same as in other MF approach pT| , ^2[ , 
but their significance is completely different being defined 
with respect to a different scaling field. 

We have obtained a complete characterization of the 
critical properties of the sandpile model. The critical 
state arises due to the fine tuning of the driving rate and 
the dissipation. This condition is enforced implicitly in 
the BTW model by imposing time scale separation and 
dissipation only through the boundaries, which makes h 
and e equal to zero in the thermodynamic limit. In this 
formalism SOC appears as a special case of non equilib- 
rium critical phenomena, with the only peculiarity that 
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the critical parameters are zero. The same MF analysis 
applied to the forest-fire model leads to similar conclu- 
sions |l(]j24|]. Contrary to sandpile models in the forest- 
fire model no conservation laws are present and MF ex- 
ponents are not correct in low dimensions. It would be 
desirable to treat extremal models jl0| in a similar way, 
emphasizing the role of the driving rate. We hope this 
will clarify the precise significance of SOC in the frame- 
work of non equilibrium critical phenomena. 

The Center for Polymer Studies is supported by NSF. 
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FIG. 1. The density of active site in the BTW model with 
border dissipation as a function of the driving rate h. 
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FIG. 2. The density of active sites for the BTW model 
with bulk dissipation e and periodic boundary conditions as 
a function of h/t for L=64. 
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